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ABSTRACT 

Aims. More than half of all active galactic nuclei show strong photoionised outflows. We use a massive multiwavelength monitoring 
campaign on the bright Seyfert 1 galaxy Mrk 509 to constrain the location of the outflow components dominating the soft X-ray band. 
Methods. Mrk 509 was monitored by XMM-Newton and other satellites in 2009. We have studied the response of the photoionised 
gas to the changes in the ionising flux produced by the central regions. We used the five discrete ionisation components A-E that we 
detected in the time-averaged spectrum taken with the RGS instrument. By using the ratio of fluxed EPIC-pn and RGS spectra, we 
were able to put tight constraints on the variability of the absorbers. Monitoring with the Swift satellite started six weeks before the 
XMM-Newton observations. This allowed us to use the history of the ionising flux and to develop a model for the time-dependent 
photoionisation in this source. 

Results. Components A and B are too weak for variability studies, but the distance for component A is known from optical imaging 
of the [O in] line to be about 3 kpc. During the five weeks of the XMM-Newton observations we found no evidence of changes in 
the three X-ray dominant ionisation components C, D, and E, despite a huge soft X-ray intensity increase of 60% in the middle of 
our campaign. This excludes high-density gas close to the black hole. Instead, using our time-dependent modelling, we find that the 
density is very low, and we derive firm lower limits to the distance of these components. For component D we find evidence for 
variability on longer time scales by comparing our spectra to archival data taken in 2000 and 2001, yielding an upper limit to the 
distance. For component E we derive an upper limit to the distance based on the argument that the thickness of the absorbing layer 
must be less than its distance to the black hole. Combining these results, at the 90% confidence level, component C has a distance 
of > 70 pc, component D is between 5-33 pc, and component E has a distance > 5 pc but smaller than 21^100 pc, depending upon 
modelling details. These results are consistent with the upper limits that we derived from the HST/COS observations of our campaign 
and point to an origin of the dominant, slow (v < 1000 km s _I ) outflow components in the NLR or torus-region of Mrk 509. 

Key words. Galaxies: active - quasars: absorption lines - X-rays: general — X-rays: galaxies — Galaxies, individual: Mrk 509 



1. Introduction 

It has been known for a long time that accretion onto a super- 
massive black hole plays a key role in the physics of active galac- 
tic nuclei (AGN). While outflows have alread y been detected us- 
ing ea rly optical observations of NGC 4151 (Anders on & Krafil 
1969), their widespread occurrence and role was only recog- 
nised after high -resolution UV spectrographs became available. 
ICrenshaw et al.1 d 19991) found that 60% of all Seyfert 1 galaxies 
show UV absorption lines. While intrinsic X-ray absorption had 



already been found in the early 1980s, including the first evi- 
dence of absorption from ionised gas (Halpern 1984), the proof 
that this gas is outflowing came fr om the first high-re solution 
X-ray spectrum of a S e yfert g alaxy dKaastra et al.ll200d) . We re- 
fer to ICrenshaw et alj d2003 ) for a more extensive overview of 
outflows. 

Measurements of absorption line spectra yield reliable in- 
formation on different aspects of the outflow such as dynamics 
(through the Doppler shifts and line profiles) and ionisation state 
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(through intensity ratios of lines from different ions). However, 
they do not provide a key ingredient, the distance r of the ab- 
sorbing material to the central black hole. The basic problem is 
that the ionisation state is described by the parameter £ = L/nr 2 
with L the ionising luminosity between 1-1000 Ryd and n the 
hydrogen density. With L and £ known from observations, only 
the product nr 2 is known. This degeneracy has led to consider- 
ing many different locations for the outflow, including the halo 
of the host galaxy (10 kpc), the inner galactic disc (0.1-1 kpc), 
the narrow emission line region (10-500 pc), and the interme- 
diate line region (0.01-10 pc ). The recently discu ssed, so-called 
ultra-fast outflows (see e.g. Tom besi et al.l l2010, and references 
therein) may be even further inwards. 

There are several methods of estimating r. The most direct 
and easiest method is to use direct im aging. A good exam ple 
are the Doppler maps of [O m] ,15007 bv lphillips etal] d 1983b in 
Mrk 509. However, this method can only be applied for the most 
nearby sources and for the most extended low-ionisation gas. 
Instrumental capabilities are insufficient to resolve more cen- 
trally concentrated gas. 

Other methods rely on obtaining the density n and deriving 
the distance r through using the ionisation parameter. The most 
direct way is to use density-sensitive lines from meta-stable lev- 
els. This method has been applied in a few cases using different 
UV lines from meta-stable levels in various sources. The pub- 
lished results display a broad range of dis tances for different 
sourc es, from <0.1 pc ( Krae mer et al.ll2006|). 25 p c ( Gabel etal] 
120051) to 3 kpc dMoe et al.l2009h . In X-ravs. kaastra et al.f l2004) 
pioneered this method theoretically using K-shell absorption 
lines from O v, but owing to the limited signal from the target 
that was studied, Mrk 279, they could not obtain firm detections 
of these lines. For X-ray binari es, however, metas table lines have 
been found in GRO J 1655^10 dMiller et alj|2008l) . but they sam- 
ple a density range (~10 20 m~ 3 ) that is too high for AGN. 

Alternatively, variations in the ionising luminosity of an 
AGN will induce changes in the ionic composition of the sur- 
rounding gas, hence changes in the transmitted spectrum. How 
fast the gas responds depends on the recombination time scale, 
which is a function of density (for absorption the light-crossing 
time is not relevant). By measuring the delay (or lack of delay) 
of an absorption component relative to the luminosity variations, 
this recombination time scale can be determined and from that 
the density n and distance r are deduced. The method has been 
applied to several sources, again yielding a broad range of dis- 
tanc es (see for reference s our introduction in paper I of this se- 
ries, |Kaaiiaitin|20nB . 

In this paper we follow the last method. In 2009 we have 
performed a large multiwavele ngth campaign on on e of the best 
targets for this study, Mrk 509 ( Ka astra et al.l201 lbh . The source 
has a typical variability time scale of a few days, with suffi- 
ciently large amplitude. The first step of our analysis was to 
characterise the time-averaged properties of the absorber using 
the combined RGS spectrum of all observations taken by XMM- 
N ewton during our cam paign. This has been described in detail 
bv lDetmers et al. | (|201 ll). Using the m ultiwavelength monitoring 
as described bv lKaastra et alj d201 lblh it is then possible to pre- 
dict the variations of the transmission of the outflow for a grid of 
possible distances r. Comparing these predictions with the ob- 
served spectra then yields distances to the outflow components. 
In our comparison, we have focussed on the data obtained with 
the pn camera of XMM-Newton, because that has the highest 
sensitivity to variations in the 1 keV band. The MOS data could 
not be used for this comparison because of a small amount of 



pile-up. However, we also obtained useful constraints from the 
RGS spectra at longer wavelengths. 

The structure of this paper is as follows. First we summarise 
the ionisation structure of the outflow, investigate the spectral 
signatures of the expected variations of the transmission, and 
introduce the data sets that we use (Sect. [2]). As a next step 
(Sect. H]i we compare the data to the predictions for the high- 
density limit (instantaneous response of the absorber to contin- 
uum variations). Afterwards we describe our method of treat- 
ing the more general case of time-dependent photoionisation 
(Sect. |5J. In Sect. [6] we compare the data from our 2009 cam- 
paign to archival data to obtain constraints from the long-term 
variability. We discuss our results briefly in Sect. [7] 

We use lcr (68% confidence) uncertainties throughout this 
paper, unless stated explicitly otherwise. 

2. Data 

2.1. Fluxed spectra 

We have made fluxed spectra for each of the ten observations as 
well as for the combined spectrum for each of the instruments 
available to us. The need to use fluxed spectra is explained in 
more detail in Sect. [3] and is related to the weakness of the ex- 
pected signals. More details a bout the observations and the ob- 
servation log are described bv lKaastra et alj d2011bl) . 

2.1.1. EPIC pn spectra 

The pn dStriider et al.ll20Q TT) count spectra were extracted as de- 
scribed by Ponti et al. (201 1, in preparation). We have produced 
fluxed spectra from these count rate spectra using the program 
epicfluxer that is available in the public SPE^Q distribution. 
Our procedure is as follows. We rebin the count spectra energy- 
dependent to about 1/3 FWHM, ignore the parts below 0.2 keV 
and above 12 keV, and fit the count spectra using a cubic spline 
(the spin model of SPEX) on a logarithmic energy grid between 
0.2-12 keV with 100 grid points. The resolution of this grid is 
about 2%, and it matches the spectral resolution of pn well. All 
fits are statistically acceptable. We then take the corresponding 
model photon spectra (also as splines) and convolve them with 
a Gaussian redistribution function. The FWHM of this Gaussian 
is energy-dependent and corresponds to the spectral resolution 
of the instrument at each energy. In this way we have a model 
spectrum at EPIC resolution that describes the data perfectly. 
Any remaining residuals that are visible in a plot of the original 
fit to the count spectra with the spin model, come mostly from 
the statistical noise of the data. To take that noise into account 
and also to correct for any remaining weak residuals caused by 
the limitations of the spline photon spectrum approximation to 
the true photon spectrum, we scale the convolved spline model 
spectra with the ratio of the observed data to the predicted data 
at each energy. 

Fig. Q] shows an example of the fluxed spectrum of the com- 
bined 2009 data. For our analysis we used the data points, but 
the plot shows the excellent agreement between the response- 
convolved spline model and the data. The spectrum clearly 
shows the Fe-K line near 6.4 keV, as well as a dip caused by 
the absorbing outflow between 0.8 and 1.0 keV. Some of the 
features between 1.8 and 4 keV may be caused by remaining 
small calibration uncertainties at a level of a few percent, and are 
possibly associated to the gold edge of the mirror or the silicon 

1 www.sron.nl/spex 
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Fig. 1. Fluxed pn spectrum for the co-added 2009 spectrum of 
Mrk 509. The spectrum has been multiplied by E 2 for clarity in 
the display. 
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Fig. 2. Ratio of the 'fluxed' pn count spectrum to the fluxed pho- 
ton spectrum of Mrk 509. 



edge of the detector. However, whether these features are instru- 
mental or have an astrophysical origin is not important for the 
present purpose, as they disappear in the ratio of different spec- 
tra. Below 0.4 keV the flux drops rapidly owing to the Galactic 
absorption. Since the spectral energy resolution in that region of 
the spectrum is low anyhow and we do not expect to see signif- 
icant variations in the absorption components in that band, it is 
not important for the present paper, so we ignore the spectrum 
below 0.5 keV in the rest of the paper. 

Furthermore, we take the following into consideration. 
Fig. |2] shows the ratio of the nominal count spectrum 
(counts s keV -1 ) divided by the nominal effective area at each 
energy, to the fluxed spectrum produced by the epicUuxer pro- 
gram as described above. In particular, at lower energies most 
of the counts correspond to redistributed photons of higher en- 
ergies. Because the intrinsic photon spectrum becomes weak at 
low energies owing to interstellar absorption, the strong rise of 
the ratio below 0.35 keV, up to ~ 4 for E = 0.2 keV, is produced. 
The negative dip around 0.45 keV is caused by loss of photons 
of 0.45 keV to counts at energies below 0.45 keV. 



In the data, an excess of 1% above the model at 0.5 keV 
would actually be an excess of 1 % counts. With the ratio of ~ 0.9 
at 0.5 keV from Fig. [2] this corresponds to an excess in flux of 
about 1.1%. Our procedure would thus give a systematic bias 
of 0.1% for the above example. However, the typical statistical 
uncertainty of the spectra in this region is much smaller than 
the 1% of the example (<0.4-0.5%). Furthermore, the statistical 
deviations for neighbouring bins are independent and can have 
different signs, so their average value over a certain band of bins 
is even lower. We conclude that these effects on our analysis can 
be ignored safely. 

The pn is known to be extremely stable over time. Here we 
quantify this. Stability can only be measured using unvarying 
sources. Such sources are relatively rare, but perhaps the best 
candidate is the isolated neutron star RXJ 1856.6-3754. This 
object is being used as a regular calibration target for XMM- 
Newton. The best band for checking stability is the 0.3-0.5 keV 
band, given the soft spectrum of the source. We use the prod- 
ucts that are available with the calibration review tool on the 
XMM-Newton web page. The 15 pn observations of this source 
in the 0.3-0.5 keV band have fluxes with a standard deviation of 
3.8% around the mean value. We are more interested in relative 
differences rather than absolute fluxes, so we have also deter- 
mined the mean and standard deviation for the fit residuals aver- 
aged over a 0.05 keV wide band around 0.3 keV. These residuals 
have a mean value of 2.4% (consistent with the nominal cali- 
bration uncertainty), but more importantly the scatter (standard 
deviation) for individual observations around this mean is only 
~ 0.5%. The true number may be even smaller because this num- 
ber agrees fully with the uncertainty due to the photon statistics 
in this band (it is limited by the exposure times of the calibration 
observations). The number implies that relative time variations 
in narrow bands can be measured down to the ~ 0.5% level or 
better, provided the source has enough photons. 



2.1.2. EPIC MOS spectra 



The MOS dTurner et al.ll200l1) data cannot be used for our pur- 
pose. In contrast to the pn data, the ratios between MOS spectra 
often show weak breaks around 0.7 keV, with deviations below 
that energy from the extrapolated power law from the 0.7-1.5 
keV band and also deviations from the ratio obtained from the 
simultaneous pn data. This is caused by a combination of small 
amounts of pile-up, small imperfections in the MOS redistribu- 
tion function at low energies, and uncertainties near the instru- 
mental oxygen edge. This softening or hardening near 0.7 keV 
makes detection of the expected narrow Gaussian differences in 
transmission harder, because at least for the ionisation compo- 
nents C and D, the centroids of the Gaussians are near this en- 
ergy (0.716 and 0.820 keV, Tabled. The signal-to-noise ratio of 
the combined MOS data in the relevant band is also about two 
times lower than for the corresponding pn data. This is partly 
due to the smaller extraction radii for the MOS spectra (20") as 
compared to the pn spectra (45"). For these reasons we do not 
include the MOS data in our subsequent analysis. 

2.1.3. RGS spectra 



The data reduction of the R GS (Iden Herder et al.l |200~TT) spec- 
tra was described in detail by Kaa stra et al.l ( 201 Tab . The spectra 
were taken in the so-called multi-pointing mode. For each ob- 
servation, one of five offset pointings in the dispersion direction 
was chosen. This has the advantages that (i) the effects of spec- 
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Fig. 3. lcr relative uncertainty on the individual ~60 ks RGS 
spectra of Mrk 509 per 0.058 A bin. 



tral bins with lower exposure from the excision of bad pixels is 
strongly reduced in the stacked RGS spectrum, and (ii) that the 
gaps between CCDs are filled with data. For individual spectra, 
this advantage is lost, although it is partly recovered by the com- 
bination of RGS1 and RGS2 data and first- and second-order 
spectra. 

We are searching for variation in the depth of absorption 
lines, so we have rebinned our spectra by a factor of 6 to 0.058 A, 
approximately the resolution (FWHM) of the RGS. In Fig. [3] 
we show our sensitivity to detecting changes in individual RGS 
spectra. The curve is derived from a 50-point spline fit to the 
median statistical uncertainties of the ten individual spectra. It 
ranges between 3 and 6% in the 10-30 A band, with some max- 
ima (7-9%) at the two ends of the band and in correspondence 
with the missing CCDs. 



Table 1. Ionisation components used in the present analysis 



component 


logf 


A'h 6 


A 


-0.33 ± 0.49 


0.23 ± 0.09 


B 


0.71 ±0.12 


0.84 ±0.10 


C 


2.01 ±0.02 


4.8 ± 0.4 


D 


2.79 ± 0.06 


5.7 ±0.9 


E 


3.60 ± 0.27 


54 ±73 



Notes. (fl) Ionisation parameter in 10 9 Wm as used throughout this 
paper. (A) Total hydrogen column density in 10 24 irT 2 . 
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Fig. 4. Transmission for the five ionisation components A-E at 
the resolution of the pn instrument (solid lines). Dotted lines: 
the same but for log£ 0.1 dex higher. Dashed lines: ratios of the 
transmission for the higher of the above ionisation states to the 
lower ionisation states. We use the same colour coding for each 
component. 



3. Description of the outflow and analysis method 

3.1. Ionisation structure of the outflow 

The time-aver a ged RGS spectrum of Mrk 509 as analysed by 
iDetmers et al.l (I201 ll) shows a number of discrete ionisation 
components, labelled A-E (see Table Q] taken from the absorp- 
tion measure distribution fit). The discrete nature of these com- 
ponents could be proved - at least for components C and D - 
by applying a continuous model for the column density versus 
ionisation parameter That continuous model shows two strong 
peaks in absorption measure at the position of components C and 
D. The statistics of the spectrum and the lack of suitable strong 
lines did not allow us to prove the discrete nature of components 
A, B, and E, but there is also no strong evidence against the dis- 
crete nature of these components, so for o ur present analysis w e 
started with the five components given bv lDetmers et al.l (I201 ll) . 
We use the parameters from his Table 6, given here in Table [1] 
which is based on a fit of the measured ionic column densities 
co-added over all velocities resulting in five ionisation compo- 
nents. 

Alternatively, IDetmers et al.l (1201 ll) also fitted the spec trum 
directly with the xabs model of SPEX (iKaastra et alJll99a) us- 
ing multiple discrete ionisation and velocity components. In that 
analysis, component C is split over two velocities (-40 and 



-270 kms' 1 ), with slightly different ^-values (see Table 5 of 
IDetmers et al.ll20i~fl) . In our discussion we come back to this is- 
sue in the context of discussing the uniqueness of our results. 

3.2. Transmission of the five components 

As a starting point in our analysis, we calculated the transmis- 
sion, convolved with the resolution of the pn detector and with 
the cosmological redshift applied, for all five individual ionisa- 
tion components. We focus here on the pn data because of its 
higher sensitivity at higher energies, but later we also consider 
the RGS spectra. Because Mrk 509 shows variations of 20-25% 
in luminosity (~0.1 dex), we also calculated the transmission for 
these components using the same total hydrogen column density 
but with log£ 0.1 dex higher (recall ^ = L/nr 2 ). For other, not 
too large variability amplitudes we can easily scale these calcu- 
lations. The results are shown in Fig. [4] This figure also shows 
the relative change for each component if the ionisation balance 
is adjusted, by plotting the ratio between the transmission for 
both cases. 

The figure shows that - at the pn resolution - all components 
show deviations from unity in the transmission, up to 6% for 
component C. However, the differences in the transmission be- 
tween the two cases are smaller than the deviations from unity, 
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Table 2. Approximation with Gaussian components of the 
changes in transmission due to an increase of 0. 1 dex in log 



Component 


E 


cr 


N 




(keV) 


(keV) 




C 


0.716 


0.041 


0.0113 


D 


0.820 


0.073 


0.0060 


E 


1.056 


0.098 


0.0102 



and they do not always peak at the energy of the dominant ab- 
sorption features. For instance, component D has the strongest 
absorption blend at an energy of 0.90 keV, but by increasing the 
ionisation parameter the main effect is to reduce the depth of the 
low-energy wing of the blend, resulting in a peak of the differ- 
ence spectrum at 0.82 keV. 

In general ions that have their peak concentration exactly at 
the ionisation parameter of one of the detected ionisation com- 
ponents, thereby making the strongest contribution to the ab- 
sorption, are not very sensitive to modest changes in the ionisa- 
tion parameter, since the second derivative of the ionic column 
density versus ionisation parameter is close to zero. Therefore, 
changes in the individual components can be spotted most eas- 
ily using ions that peak near the relevant ionisation parameter 
(to yield significant column densities and absorption features), 
but not exactly at this ionisation parameter (and therefore have a 
strong gradient of the derivative of the column density versus £). 

We found that the prominent peaks due to a change AT in the 
transmission for components C, D, and E are represented well by 
Gaussians: 

AT = Ne -(E-E Q ) 2 /2cr\ (1) 

The parameters for these Gaussians (for an increase of 0.1 in 
log£) are given in Table [2] 

We discuss the different components in more detail below. 
Although they affect the transmission below 1 keV at the 1— 
2% level, components A and B are expected to show changes 
in the transmission (for an increase in log£ by 0.1) no larger 
than 0.1%, so cannot be used for the present purpose of con- 
straining the spectral variations of the absorption components. 
This lack of sensitivity is due to the lower column densities for 
these low-ionisation components, combined with the poor spec- 
tral resolution of EPIC at low energies, causing relatively strong 
broadening of spectral lines. 

For this reason we concentrated on components C, D, and E. 
In their expected response to changes in the ionisation parameter, 
they show significant features with amplitudes of about a percent 
at energies of 0.716, 0.820, and 1.056 keV, respectively. Using 
our spectral model, we estimated the equivalent widths of the 
lines from all ions that contribute to these features. The results 
are shown in Table [3] 

The expected changes in component C are dominated by the 
declining column densities of Fe ix-Fe xn as a function of ioni- 
sation parameter, and also by the Rydberg series of O vn starting 
from n — 3. The ls-2p line of this ion is well separated from 
the presently considered blend. At the far blue side of the pro- 
file, the increase in Fe xvn column density enhances the contrast 
for these column densities of Fe ix-Fe xn slightly. The dominant 
lines from Fe xv-Fe xvn are separated reasonably well in energy 
from the more lowly ionised ions, and thus do not contribute 
much to the central regions of the blend. As a result, mainly the 
changes in Fe ix-Fe xiv and O vn and O vra for this component 
are relevant. 



Table 3. Strongest contributions to the dominant absorption 
troughs for components C, D, and E. 



Comp 


ion & dominant 


A" 


E» 


W, c 


W," 


AW 




transition / blend 


(A) 


(keV) 


(eV) 


(eV) 


(eV) 


C 


O vn ls-3p. . .7p 


18.63 


0.663 


0.99 


0.75 


-0.24 


C 


Feix 


16.54 


0.721 


0.17 


0.06 


-0.11 


c 


Fex 


16.36 


0.729 


0.52 


0.23 


-0.29 


c 


Fexi 


16.18 


0.738 


0.84 


0.43 


-0.41 


c 


Fexn 


16.03 


0.745 


0.69 


0.53 


-0.16 


c 


Fe xm 


15.89 


0.756 


0.60 


0.55 


-0.05 


c 


Fexiv 


15.63 


0.767 


0.27 


0.32 


+0.05 


c 


Ovm ls-3p,4p,5p 


16.01 


0.768 


1.58 


1.50 


-0.06 


c 


Fexv 


15.36 


0.777 


0.16 


0.23 


+0.07 


c 


Fexvi 


15.26 


0.786 


0.13 


0.22 


+0.09 


c 


Fe xvn 2p-3s,3d 


15.01 


0.789 


0.57 


1.07 


+0.50 


D 


Ovm ls-3p,4p,5p 


16.01 


0.766 


0.41 


0.26 


-0.15 


D 


Fexvn 


15.01 


0.794 


0.47 


0.25 


-0.22 


D 


Fe xvm 


14.20 


0.836 


1.27 


0.89 


-0.38 


D 


Fexix 


13.52 


0.881 


1.54 


1.41 


-0.13 


D 


Neix ls-2p 


13.45 


0.891 


0.20 


0.11 


-0.09 


D 


Fe xx 2p-3d 


12.84 


0.929 


0.65 


0.86 


+0.21 


E 


Fe xx 2p-3d 


12.84 


0.935 


0.07 


0.00 


-0.07 


E 


Nex ls-2p 


12.13 


0.987 


0.66 


0.46 


-0.16 


E 


Fe xxi 2p-3d 


12.29 


0.987 


0.40 


0.07 


-0.33 


E 


Fe xxn 2s-3p 


11.71 


1.038 


0.84 


0.19 


-0.65 


E 


Fe xxm 2s-3p 


10.98 


1.092 


1.01 


0.47 


-0.54 


E 


Fe xxiv 2s-3p 


10.62 


1.127 


1.06 


0.76 


-0.30 



Notes. (a) Laboratory wavelength of the strongest individual spectral 
line of the ion that contributes to the blend. (A) Predicted energy of 
the line centroid of all lines from the given ion that contribute to the 
blend, in the observer's frame (thus accounting for cosmological red- 
shift). (c) Total equivalent width of all the lines of the ion contribut- 
ing to the blend, for the parameters of the time-averaged spectrum (see 
TablefT). (</) As above, but for log if higher by 0.1 dex. (e) Difference in 
equivalent width for the case with a higher ionisation parameter com- 
pared with the time-averaged spectrum. 



Below 0.715 keV a reduction in the optical depth r of the 
Ovu edge is also visible (Fig. |4j». It decreases from 0.0112 to 
0.0079, causing the transmission above the edge to increase by 
0.33% relative to the lower energies. 

For component D, the Ovm ls-2p line, although it is the 
strongest line in this energy band, does not contribute signifi- 
cantly to the main component of the variation, because it is 150 
eV off the centroid and because the optical depth is large (2.2), 
causing moderate saturation effects, hence less sensitivity to col- 
umn density variations. We have included Fe xx for component 
D in Table [3] because this ion is expected to increase its col- 
umn density as a response to an increase in ionisation parameter, 
resulting in a shift of the centroid of the varying transmission 
component to slightly lower energies. 

The changes in component E are dominated by decreasing 
column densities of Fe xxi-Fe xxiv as a function of ionisation 
parameter. Component E is also predicted to have a strong ab- 
sorption feature in the Fe-K band, mainly due to Fe xxv. At EPIC 
resolution, this line has a maximal depth of about 5 %. Because 
it is possibly blended by emission features of the same ion, it is 
not possible to measure its equivalent width. If component E re- 
sponds to the continuum changes, then the line's peak intensity 
will change by about 1 %, which is too small to detect since the 
signal-to-noise ratio per observation and per resolution element 
(FWHM) is only 50 at those high energies. Therefore we do not 
consider the Fe-K band further here. 
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For component E, there is a significant uncertainty of 0.27 in 
th e value of log £, for the time-averaged spectrum as determined 
bv lDetmers et al.l d201 lb . This is because of a strong correlation 
between the total hydrogen column density and The conse- 
quences of this on our analysis are discussed in more detail in 
Sect. 13 For components C and D, the uncertainties are much 
smaller. However, because we use in this paper changes rela- 
tive to the time-averaged spectrum, uncertainties in the absolute 
values of £ and Nh for the time-averaged spectrum are less im- 
portant. 




3.3. Measuring transmission changes of the outflow 

The continuum flux during the XMM-Newton campaign of 
Mrk 509 varied significantly at a level of the o rder of 25%, 
corres ponding to 0.1 dex, in the hard X-ray band (Kaastra et al. 
1201 lbl) . In the soft X-ray band we even see enhancements up to 
~60% (see Fig. [T2l later in this paper). If any of the ionisation 
components of the outflow responds rapidly to the changing ion- 
ising continuum, then we expect to see the changes in the spec- 
trum as displayed in Fig. |4] If on the other hand the response 
time is slow, no changes in the transmission may occur. Fig. [4] 
shows that the expected peaks for components C, D, and E are 
well separated in energy, implying that the (lack of) variations 
can be detected for each of these components individually. 

The amplitude of the expected signal is of the order of 1%. 
This is smaller than the typical systematic uncertainty in the 
EPIC effective area, which is a few percent or even slightly 
larger, and varies with energy. As a result small differences in 
the model photon spectra cannot be distinguished from imper- 
fections in the adopted effective area. Therefore, global spectral 
fitting of individual spectra at this level of accuracy is not pos- 
sible. However, the pn detector of XMM-Newton is extremely 
stable (Sect. I2.1.U . and therefore relative changes in the spec- 
trum are much easier to detect. 

We proceed as follows. We produce fluxed spectra 5, for all 
ten individual observations (z = 1,...10) as described in the 
next subsection. The ratios ry j0 b S = S ,/S j of these fluxed spectra 
are then independent of any effective area uncertainties. Also, in 
this way the effects of interstellar absorption are taken out. The 
expected ratio is r,jp re d = CiTi/CjTj, where C, is the intrinsic 
continuum for spectrum i and Tj the transmission of the absorb- 
ing component(s). If there is no change of the absorber, ry jPre d 
is just the ratio of the model continuum spectra. In order to take 
out the time variability of the continuum, we use instead 

R i j,6b, = (Si/Sj)/(C i /Cj)-l, (2) 

and compare this with 

Rij.pied = Tj/Tj - 1. (3) 

If any of the absorbing components varies over time, /?,y ;Pre d and 
the corresponding /?,j, b s will differ from zero in narrow energy 
bands, cf. our findings in Sect. 13.21 

4. Searching for variations in the outflow in the high 
density limit 

In the previous section we described the structure of the outflow, 
indicated the characteristic predicted variations as a response to 
changes in the ionising radiation, and presented the data products 
that we use to measure these variations. 

In our analysis we use three levels of sophistication. At the 
most simple level we assume that the ionising continuum varies 



0.8 - 



■ 1 —3 5 - 7 

2 4 6 

0.6 I 

10- 3 0.01 0.1 1 10 100 

Energy (keV) 

Fig. 5. Ratio of the ten model spectra to the time-avera ged model 
spectrum, based on the fits by Mehdino ur et al.l (I201 ll) . Only the 
underlying continua are taken into account, not the absorption 
components. The spectra as displayed are extrapolated outside 
the XMM-Newton range for illustrative purposes. 



only in flux, but not in shape. In the case of high density, this 
causes an immediate response to the flux variations with the 
characteristic signatures as presented in Fig. [4] For gas at a fixed 
density and distance to the ionising source, a change of L is then 
directly proportional to a change in 

However, t he above situation is not realistic, 

iMehd ipour et al] (|201 lb have shown that there are signifi- 
cant changes of the spectral shape during our campaign (Fig. [5] 
Sect. I4.1l i. Therefore in this section we study a model where 
the gas has a high density and hence responds immediately 
to variations in the ionising flux, but we take the variations 
in spectral shape as measured simultaneously with the EPIC 
and OM instruments of XMM-Newton into account. This is an 
important limit for the more general model of time-dependent 
ionisation at lower densities that we discuss in Sect. [5] 

4. 1 . Ionising continuum 

We follow Mehdin our et al.l (I201 lh . who model the continuum 
as the sum of four components: a disc blackbody (dbb), a soft 
excess parametrised by a Comptonised disc component (se), a 
power law (pi) and an Fe-K line (fe, see their Table 5). The iron 
line does not contrib ute much to the ionising flux, but for consis- 
tency with the fits by Mehdip our et al.l (I201 ll) we include it in our 
model. However we ignore the wea ker narrow and broad X-ray 
emission lines that are modelled bv lDetmers et al.l d201 lh in the 
soft X-ray band. Those lines also have a negligible contribution 
to the ionising flux. More details about the model components 
are given in Sect. I5.21 where we discuss the time-dependent pa- 
rameters of the components. 

For each of the ten spectra, we have made runs with Cloudy 
dFerland et al.lll998l) version C08.00 with lLodders etalj (2009) 
abundances for a grid of values to calculate the equilibrium ion 
concentrations. The auxiliary program xabsinput in the public 
SPEX distribution then converts this into an input file for the 
xabs model of SPEX. The xabs model w as used in the a nalysi s 
of the time-averaged RGS spectrum by iDetmers et alJ (l201ll) . 
For each spectrum we also calculate the ionising luminosity L in 
the 1-1000 Ryd band (Tableffl). 
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Fig. 6. Data points with error bars: observed ratios of the ten spectra (labelled 1-10) to the average spectrum in the representation 
of Eq. (fj) (essentially the spectral ratios with the respective continua divided out and unity-subtracted). Solid green lines: predicted 
ratio in the representation of Eq. © for ionisation component C; dashed blue lines: same for component D; dotted red lines: same 
for component E; thick magenta line: components A-E combined. The predictions for components A and B are close to zero and 
not shown here separately. 
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Table 4. The 10 log of the ionising (1-1000 Ryd) flux for the ten 
individual spectra 1-10 divided by the corresponding flux for the 
time-averag ed spectrum, for each o f the spectral components as 
modelled by Mehdinour et all d201ll) . 



Obs 


Fe" 


PL 6 


SE C 


DBB" 


Total 


1 


-0.06 


-0.08 


-0.07 


-0.21 


-0.08 


2 


-0.05 


-0.05 


-0.06 


-0.12 


-0.06 


3 


+0.02 


+0.04 


+0.01 


-0.06 


+0.02 


4 


-0.03 


-0.07 


+0.07 


+0.18 


+0.03 


5 


+0.02 


+0.03 


+0.12 


+0.28 


+0.09 


6 


+0.01 


+0.04 


+0.02 


+0.04 


+0.02 


7 


+0.01 


+0.03 


-0.01 


-0.03 


+0.01 


8 


+0.03 


+0.05 


-0.05 


-0.12 


-0.02 


9 


+0.02 


-0.03 


+0.00 


-0.09 


-0.01 


10 


+0.01 


+0.01 


-0.04 


-0.15 


-0.03 



Notes. (fl) Iron line <fc) Power law w Soft excess <rf) Disc blackbody 



4.2. Predicted and observed changes in transmission: pn 

As a next step, we calculated for each of the five ionisation com- 
ponents A-E the transmission using the above appropriate SEDs. 
We also increased or decreased £ proportionally to the change of 
the total L from the last column of Table |4] The ratios of the ten 
transmissions to the time-averaged transmission are then repre- 
sented using Eq. [3] and compared to the observed pn data using 
Eq.|2] We show this comparison in Fig. [6] 

Observations 4 and 5 show some broad-band residuals in 
the form of a positive trend between ~ 0.8 keV and 1.5 keV, 
while observation 8 shows the opposite trend. A similar bend- 
ing is seen for these spectra in the plot of the ratio of the model 
spectra (Fig. [5}. These three spectra have their largest curvature 
just around 1 keV. Thus, the systematic trends in spectra 4, 5 and 
8 can be attributed to imperfect continuum modelling. Note that 
their amplitude of ~1% is well below the systematic uncertainty 
of the pn effective area calibration. But of course also the con- 
tinuum models that we use may be too simple at this level of 
accuracy. 

To corroborate this result, we have inspected the fit resid- 
uals of the broad-ban d fits to the ten spectra made by 
iMehdipour et alj d2011l) . Averaged over the 1.0-1.5 keV band, 
the fit residuals of the pn for spectra 4 and 5 are ~0.8% higher 
than for spectra 2, 9 and 10. Similarly, for observation 8 the fit 
residuals are ~0.4% smaller than for spectra 2, 9 and 10. Within 

0. 1%, these differences agree with what we show in Fig. [6] Thus, 
the large-scale residuals are caused by small imperfections in the 
continuum modelling of the individual spectra. Recall that the 
ratios shown in Fig. [6] are essentially the ratios of the observed 
or predicted photon spectra, but normalised to t he correspo nding 
ratio o f the pure continuum models obtained by Mehdipour et al.l 
d201 lb . In contrast to these broad-band trends, the signatures of a 
varying absorber should be narrow dips or bumps at specific en- 
ergies, which differ for each of the absorption components (cf. 
Fig. ED. 

The predicted difference in transmission for components A 
and B is less than +0.02 and +0.2%, respectively, and because 
we are not sensitive enough to detect such weak signals, we will 
not consider components A and B further. We recall that the sen- 
sitivity to detecting changes in narrow energy bands with pn is 
-0.5% or slightly better, cf. Sect. I2XT1 

Component C has the strongest predicted signals. Spectrum 

1, and to some lesser extent also spectrum 2, should have a dip 



Table 5. Comparison of the 10 spectra with models for the 
change in transmission 







Original data 




Corrected data 




Obs 


x 2a 


Ax 2b 


A* 2 * AX 2 " 


x 2a 










(C) 


(D) 


(E) 




(C) 


(D) 


(E) 


1 


60.9 


2.5 


-3.1 


2.5 


56.8 


12.2 


-1.2 


3.8 


2 


42.3 


8.4 


-0.1 


-0.6 


39.8 


7.5 


0.0 


0.0 


3 


35.9 


-3.5 


-1.2 


-0.5 


32.0 


4.6 


0.3 


0.4 


4 


59.1 


1.8 


1.4 


17.8 


29.5 


1.3 


0.0 


10.5 


5 


80.4 


19.4 


-0.7 


1.0 


33.9 


21.6 


0.7 


0.4 


6 


37.3 


1.7 


-0.8 


-0.3 


34.3 


3.9 


-0.2 


0.4 


7 


65.7 


-0.8 


-1.0 


-0.4 


63.6 


1.5 


-0.5 


0.0 


8 


58.9 


-1.1 


1.0 


7.2 


31.9 


-0.9 


0.6 


4.6 


9 


39.8 


-0.8 


-0.1 


-1.0 


31.3 


-0.2 


0.1 


-0.2 


10 


27.9 


-2.2 


-0.1 


2.8 


21.4 


-0.9 


0.1 


2.1 


S 


508.0 


25.4 


-4.8 


28.6 


374.5 


50.4 


-0.2 


22.0 







Notes. x 1 °f the zero model (A) x 2 °f me prediction for the component 
indicated in brackets minus x 2 of the zero model 



with an amplitude of — 1% at 0.73 keV, while spectrum 5 should 
have a positive excess with similar amplitude. None of this is vis- 
ible in the data. However, the predicted dip for spectrum 1 might 
be visible in spectrum 5, while the predicted bump for spectrum 
5 might be visible in spectrum 8. That would correspond to a 
delay of about two weeks. The possibility of a delayed response, 
corresponding to lower density and hence more distant gas will 
be elaborated further in Sect. [5] but we conclude here that the 
pn data are clearly inconsistent with an immediate response of 
component C to continuum variations. 

The predicted changes for component D are too small to be 
detectable with pn. This is somewhat surprising because from 
Fig. [4] one might conclude that the predicted variations should 
have about half the amplitude of the variations in component C. 
However, in Fig.|4]it was assumed that the luminosity variations 
are achromatic, while in the present model we take the more re- 
alistic continuum shape variations into account. Because the soft 
X-ray flux varies more strongly than the hard X-ray flux, the pre- 
dicted amplitudes for the more highly ionised components D and 
E are smaller than that of component C. The strongest predicted 
dip, -0.35% around 0.82 keV for spectrum 1, is comparable to 
the statistical uncertainty of 0.3% in that band, hence not de- 
tectable. 

For component E, our model predicts a dip of -1.6% for 
spectrum 4 at 1 .06 keV, in contrast to a positive peak at the same 
energy of — h0.5% for spectra 6, 8, and 10. There is not much 
evidence of the dip in spectrum 4, although some care should be 
taken because of the small systematic trends noticed earlier for 
this spectrum, related to the strong curvature of the underlying 
continuum. 

4.2.1 . Nominal significance of the variations 

The significance of our findings is somewhat difficult to assess 
quantitatively, due to the small systematic continuum deviations 
mentioned earlier. We have attempted to quantify this in the fol- 
lowing way (Table|5]). We calculated the nominal^- 2 of the data 
i?obs shown in Fig. [6] compared to two models. The first model 
(labelled "zero model") in the table compares the data to the 
model where the change in transmission is zero (R pK d = in 
terms of Eq. [3}. This corresponds formally to the limit for the 
gas density n — » 0. We also calculated ;f 2 for the models with in- 
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Fig. 7. Example of continuum adjustment for spectra 4 (top 
panel) and 5 (bottom panel). Data points with error bars: ob- 
served ratio of spectra 4 and 5 to the average spectrum in the 
representation of Eq. (0, same as in Fig. [6] Solid black line: cu- 
bic polynomial fit to R over the 0.4-3.2 keV band as described 
in the text. Solid green line: predicted R in the representation 
of Eq. © for component C, but added here to the continuum 
approximation given by the black line; dashed blue line: same 
for component D; dotted red line: same for component E; thick 
magenta line: components A-E combined. 



stantaneous response for components C, D, and E (the green, 
blue and red curves in Fig. |6). This corresponds formally to 
n — > oo. For comparison with the zero model we subtract x 2 
for that model, thus showing the difference in x 2 - We do this 
for each individual observation. The table clearly shows that the 
model with no response at all is better than the models with an 
instantaneous response of the absorption components. The num- 
bers agree with what we described qualitatively before. While 
the sensitivity to detecting changes in component D is low, an 
instantaneous response for component C would make the total 
X 2 25 A higher than the zero model. An almost similar conclu- 
sion can be drawn for component E. 

To test how these results depend on the systematic uncer- 
tainties in the continuum, we re-computed these x 2 values using 
corrected data (the rightmost four columns of Table[5]l. The cor- 
rection consists of a cubic polynomial in log E that we fitted to 
the data between 0.4-3.2 keV, thus with a typical resolution of a 



16.5 
Wavelength (A) 

Fig. 8. Scaled flux ratios compared to the time-average spectrum 
(eqs. [2|-[3]l for RGS. Upper panel: data and model for spectrum 
5; Lower panel: same for the combined spectra 1 & 2. The data 
have been rebinned by another factor of two to 0. 1 16 A for clar- 
ity of display. The models are for component C only. For com- 
parison with other figures, the band that is plotted corresponds 
to 0.69-0.83 keV. 



factor of 2 in energy. An example of such a fit for spectrum 5 is 
shown as the black line in Fig. [7] The correction is not unique, 
and there is some danger that it washes out a small part of any 
true variations of the absorber. Keeping this in mind, the correc- 
tion lowers the overall values of x 2 compared to the uncorrected 
case as expected. It does not alter our conclusion about poorer 
fits with models with an instantaneous response compared to 
models with no response. Thus, we confirm our more qualita- 
tive conclusions given before. 

4.3. Predicted and observed changes in transmission: RGS 

At energies below 1 keV, RGS becomes more sensitive to the de- 
tection of weak lines than pn, thanks to its higher spectral resolu- 
tion that compensates for the lower effective area. Therefore we 
also consider the predicted instantaneous response to continuum 
variations for the absorption components as measured by RGS. 
We consider only features where the RGS has enough sensitivity 
to detect them (Fig. [3}. 

As for the pn, variations in components A and B are below 
the detectability limit so are not considered further. For com- 
ponent C, the 15-18 A band is most promising for detecting 
changes (Fig. |H). Our model predicts the strongest contrast be- 
tween spectra 1 and 2 on one side, which have almost the same 
transmission, and spectrum 5 that corresponds to the peak of the 
soft X-ray emission. 

Two prominent features contribute to the difference shown 
in Fig. [8] First, the model predicts a strong change in a broad 
band between 16.5-17.5 A. This is the same blend as the pre- 
dicted blend in the pn data around 0.73 keV (Fig. |6}. The main 
contributions are from to M-shell transitions of Fe ix-Fe xm, cf. 
Table[3] and the O vn edge. Averaged over the 16.5-17.5 A band, 
the model predicts an increase of R by +0.043 from observa- 
tions 1 & 2 to observation 5. The data show a decrease instead 
of -0.028 + 0.011. Thus, at the 6.5cr confidence limit we can 
exclude an instantaneous response for component C. 
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The second feature is the Fe xvn resonance absorption line 
(rest-wavelength 15.01 A, observed at 15.5 A due to the cos- 
mological redshift). Averaged over two adjacent bins, the model 
for component C predicts a decrease of -0.060 from observa- 
tions 1 & 2 to observation 5, while the data shows a slight in- 
crease of +0.05 + 0.03. A part of this increase can be explained 
by contamination by component D, which also produces Fe xvn. 
Component D alone produces an increase of +0.024. Because 
the contributions from components C and D to the total Fe xvn 
column density are 59 and 41%, respectively (cf. Table 4 in 
iDetmers et al.ll20i"il) . the combined predicted net decrease in the 
Fexvn line is -0.026. This is off by 2.5cr from the observed 
value. Thus, the Fe xvn line is also not in good agreement with 
an instantaneous response. 

For component D the largest changes are expected between 
spectrum 1 and 5, for the Ovm Lye line at 19.6 A. Averaged 
over two bins, the observed increase of 0.02 + 0.05 agrees with 
the prediction of +0.04, but of course also with the hypothesis 
of no changes. 

Finally, for component E there is a group of features in the 
1 1.3-12.8 A band where the model predicts a significant change, 
from a value of -0.015 for spectrum 4 to an average of +0.007 
for observations 6, 8, and 10. This increase by +0.022 is signif- 
icantly higher (at the 2.5<x level) than the observed difference of 
-0.003 ± 0.010. Alternatively, fitting a linear function to the ob- 
served /^-values for all ten observations to the predicted values 
yields a slope of +0.18 + 0.40, 2<x off from the expected value 
of 1 for an immediate response. This strengthens our conclusion 
from the pn data that also component E does not show much 
evidence of an instantaneous response. 

In the next section we investigate whether the observations 
agree with a lack of response or with a delayed response. For 
this, we need to solve the time-dependent ionisation balance 
equations. 

5. Time-dependent photoionisation 

5.1. Introduction 

The time evolution of the ion concentrations in the outflow de- 
pends on the gas density. For high density, the gas responds fast 
to changes in the ionising flux and restores ionisation equilib- 
rium at the new flux level. For low density, collisions between 
free electrons and ions are rare, and because of that it takes a 
long time to restore equilibrium. When the photoionising flux 
varies faster than this response time, the plasma may be continu- 
ously out of equilibrium and it is at best in equilibrium with the 
long-term average flux ; on time scales longe r than the response 
time of the plasma. See Kroli k & Krissldl995l Sect. 6) for a more 
extensive discussion of these effects. 

For strong enhancements of the flux, the i onisation may 
proce ed faster than recombination (see Fig. 5 of iNicastro et al.l 
1999, for a good example of this), but when the flux does not 
vary by orders of magnitude, as in the case of Mrk 509, the ion- 
isation and recombination time scales are of the same order of 
magnitude. 

In this section we present how w e model the tim e -depe ndent 
ionisation following the method of INicastro et al.l (U999). We 
first determine the history of the ionising spectral energy distri- 
bution (SED) for our observations. Starting with the first Swift 
observation, for which we assume ionisation equilibrium, we 
calculate for each outflow component the time evolution of the 
ion concentrations, for a grid of values for the gas density. We 



then calculate the corresponding X-ray transmission for the dif- 
ferent components, and compare those to the observed spectra. 
The density for each component (or a lower- or upper limit) is 
then determined from the best match of a model spectrum to the 
observed spectrum. 

5.2. Time-dependent SEDs 

As in Sect. 14. 11 we use the detailed spectral fits of the ten XMM- 
Newton observations from Mehdi pour et al.l d201 ll) . They model 
the continuum with four components: 

1. disc blackbody (DBB, SPEX model dbb), representing the 
direct emission from the accretion disc. Free parameters are 
the emitting area Adbb and the temperature Tdbb\ 

2. A soft X-ray excess (SE, SPEX model comt), representing 
Comptonised emission from the disc in a warm plasma. Free 
parameters are the seed temperature of the photons (taken to 
be equal to Tdbb), the plasma temperature T e , optical depth 
r, and normalisation iV comt ; 

3. power law (PL, SPEX model pow), with free parameters the 
photon index F and the normalisation A^l; 

4. Fe-K line (Fe, SPEX model gaus), approximated by a 
Gaussian and characterised by its energy Ep e , normalisation 
A^e, and Gaussian width <x. 

We apply a high-energy cut-off to the s pectrum, based on the 
INTEGRAL data modelling, as shown in lKaastra et alJ d2011bl 
Fig. 3). The cut-off is represented by a multiplicative factor to 
the power-law component, and the shape of the cut-off is kept 
constant in time, because more details are lacking. 

The bes t-fit parameters for the D BB, SE, and PL components 
are given in Mehdip our et al.l (1201 ll Table 6). The parameters for 
the iron line are not given in that paper, but the effects of the iron 
line on the present analysis are weak, because its ionising flux is 
only ~ 0.5 and 0.2% of the ionising flux of the power-law com- 
ponent and the total ionising flux, respectively. For simplicity, 
we assume that its flux is constant at the average value for the 
ten XMM-Newton observations, 0.76 phrrT 2 s _1 , corresponding 
to 2.07 x 10 35 W in the source rest frame. Our continuum model 
is therefore described by seven free parameters that vary as a 
function of time. 

Because the ionisation state of the gas also depends on 
the history of the ionising flux before the start of the XMM- 
Newton observations, we also need to characterise that part of 
the light curve. For that purpose we have ob tained Swift obser- 
vations as described by lKaastra et al.l (l2011bl) and elaborated by 
iMehdipour et alJ d2011l) . While for the ten XMM-Newton obser- 
vations we had high signal-to-noise pn spectra and OM observa- 
tions in all filters to constrain the SED, for the Swift observa- 
tions, which had short exposure times, we have in most cases 
only UV fluxes for the UVM2 filter and broadband X-ray fluxes 
at 0.3 and 4 ke V with a typical uncertai nty of 0.05-0. 10 dex (see 
Figs. 2 and 8 of [Mehdipo ur et alj|2011h . 

These three fluxes (UV, 0.3 and 4 keV) per Swift observa- 
tion are in principle insufficient for deriving the seven continuum 
parameters for the SED. However, there are tight correlations 
between these parameters, and therefore we can make a good 
proxy for the SED based on the UVM2 and the 4 keV fluxes 
alone, as we show below. We use the higher signal-to-noise and 
more detailed XMM-Newton observations to derive these cor- 
relations. For the Swift observations without the UVM2 filter, 
we estimate the UV flu x from the other filters as described by 
IMehdipour etalJd201lh . 
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UV flux at 2310 A(10- 14 Wnr 2 ) 

Fig. 9. Correlation between the photon index of the PL compo- 
nent with the UV flux for the ten XMM-Newton observations. 
The corners of the diamonds correspond to the nominal statisti- 
cal uncertainty on both quantities. 



For the power-law component, the photon index F correlates 
remarkably well with the UV flux (Fig.|9j- The PL is then com- 
pletely determined from the UV flux and the 4 keV X-ray flux, 
which has no contribution from the other softer continuum com- 
ponents. 

For th e soft excess compon ent, the plasma temperature as 
given by Meh dipour et alj d201 ll) shows little variation, and thus 
we keep it constant at 0.20 keV. The normalisation N coml and 
optical depth r of this component can be approximated as a lin- 
ear function of the UV flux (Fig. [TOb . and as stated before, the 
temperature of the seed photons is tied to the disc blackbody. 

Finally, for the disc blackbody component, there is a correla- 
tion between its temperature and the UV flux (Fig.fTTIi. We kept 
its normalisation constant to 3.34 x 10 25 m 2 , the mean value for 
the ten observation s, and within the error bars, consistent with 
the reported values (Mehdipou r et alfcoi ll Table 6). 

Using these correlations and the UVM2 and 4 keV flux, we 
constructed the time-dependent SEDs for our full 100 day cam- 
paign. Obviously, for the XMM-Newton obse rvations we took 
the SEDs directly from lMehdipour et al.l (1201 ll) and we reported 
them in Sect. 14.11 A comparison of the ionising fluxes deter- 
mined directly from the XMM-Newton data with the ionising 
fluxes determined through the above correlations shows a scatter 
of- 10%. 

We show the ionising fluxes of the main components (PL and 
SE) in Fig. [12] It is seen that the luminosity before and after the 
XMM-Newton observations is lower than during these XMM- 
Newton observations. Most of the variability is caused by the 
soft excess component, although the power-law component also 
varies. The ionising luminosities of the disc-blackbody and iron 
line are much lower than the luminosities of the PL and SE com- 
ponents. 

5.3. Solution of the time-dependent rates 

The time dependence of the relative density n, of an ion i of a 
given species is given by 

dtij 

— = -riin e a^i + n i+ in e ai - n-,!-, + n,-i/,-i. (4) 




16 17 18 19 20 

UV flux at 231 A (1 0- 14 Wrrr 2 ) 




16 17 18 19 

UV flux at 2310 A (10- 14 Wnr 2 ) 



Fig. 10. Correlation between the normalisation N comt and opti- 
cal depth t of the soft excess (SE) component with the UV flux 
for the ten XMM-Newton observations. The corners of the dia- 
monds correspond to the nominal statistical uncertainty on both 
quantities. 




16 17 18 19 

UV flux at 231 A (1 0- 14 Wnr 2 ) 



Fig. 11. Correlation between the temperature of the DBB com- 
ponent and the UV flux for the ten XMM-Newton observa- 
tions. The corners of the diamonds correspond to the nominal 
statistical uncertainty on both quantities. 
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Fig. 12. Ionising luminosity (1-1000 Ryd) for the power-law 
component (blue), Comptonised component (red), and total 
spectrum (black). Filled triangles: based on Swift data; open cir- 
cles: XMM-Newton observations. 



Here n e is the electron density, and a-, is the recombination rate 
from stage i + 1 to stage i and is in general a function of the tem- 
perature T. The ionisation rate from stage i to i + 1 is indicated 
by Obviously, in Eq. [4] ionisations from bare nuclei and re- 
combinations onto neutral atoms should be excluded as they are 
non-existent. 

Equation (0]i forms a set of coupled ordinary differential 
equations. We solve these numerically by using a Runge-Kutta 
method with adap tive stepsiz e control, subroutine odeint from 
iPress et all (119921) . Following iNicastro et all dl999l) . we adopt a 
constant electron density n e as a function of time. 

The recombination and ionisation rates are interpolated lin- 
early between the values for the 29 different Swift and XMM- 
Newton observations. We also tried quadratic or spline interpo- 
lation, but that sometimes leads to unwanted overshooting, in 
particular near the start and the end of the XMM-observations 
where a few Swift observations are close in time to those by 
XMM-Newton. Moreover, such higher order schemes would 
hide the intrinsic uncertainties caused by the lack of data be- 
tween the 29 observations. 

The recombination and ionisation rates for the 29 epochs 
were obtained as follows. For each epoch and ionisation com- 
ponent, we calculated the e quilibrium ion conce ntrations using 
Cloud y dFerland et al.ll 19981) version C08.00, with lLodders et alj 
(2009) abundances. These equilibrium concentrations were cal- 
culated using the appropriate SED (see previous subsection). 
Furthermore, the ionisation parameters £ were scaled down rel- 
ative to t he values obtained from the average XMM-Newton 
spectrum dPetmers et al.1 [20TlT> . proportionally to the ionising 
1-1000 Ryd luminosity L. For a prompt response of the gas to 
changes in flux or SED shape, this is the appropriate scaling, as- 
suming that the gas density n and distance to the ionising source 
r remain constant in time (£ = L/nr 2 ). Using this same Cloudy 
run, we used the "punch ionization rates" option of Cloudy to ob- 
tain the corresponding ionisation and recombination rates. Since 
Cloudy takes multiple ionisation after inner-shell ionisation into 
account, the total ionisation rate obtained from the punch option 
does not return an equilibrium situation when inserted into our 
simplified Eq. |@}, which ignores multiple ionisations. To correct 
for this we have scaled the ionisation rates 7, in our code so that 




time (days) 

Fig. 13. Time evolution of the C iv concentration of component 
A, relative to the total carbon concentration. Stars: a model with 
instantaneous response to continuum variations. Coloured lines: 
predictions for absorbing material at different distances as in- 
dicated by the number near the top of the plot; these numbers 
represent the log of the distance in m. The XMM-Newton obser- 
vations were taken between the two vertical dotted lines. 



equilibrium is forced (dn,/df = in Eq. |4ji using the concentra- 
tions Hi provided by Cloudy. 

We did these calculations for a grid of distances r to the ion- 
ising source and for each of the outflow components separately. 
We assumed equilibrium for the first Swift spectrum as our ini- 
tial condition. 



5.4. Results 

5.4.1. An illustrative example 

We start here with C iv, an ion that mainly produces UV absorp- 
tion but displays a number of interesting features in the time- 
dependent ion concentration. Fig.[13]shows the predicted evolu- 
tion of the ion concentration as a function of time for the domi- 
nant C iv contributor, component A. 

For an instantaneous response, the relative C iv concentration 
varies between 0.23 and 0.34 in accordance with the continuum 
variability. For distances up to 10 17,75 m, there is not much differ- 
ence with the instantaneous response. For distances of 10 18 and 
10 18 25 m, the response becomes delayed by about three and ten 
days, respectively. These delays are related to the lower densities 
and therefore longer recombination time scales corresponding to 
the longer distances (for constant L/tj = nr 2 ). 

For greater distances, for instance 10 18 50 m, the delay be- 
comes a month but also the amplitude of the variations starts to 
decrease rapidly. The plasma is no longer able to respond to the 
fastest fluctuations of the ionising continuum. For even greater 
distances (> 10 18 75 m), the delay becomes longer than the dura- 
tion of six week of the Swift monitoring before the start of the 
XMM-Newton observations, and the amplitude becomes signif- 
icantly smaller. Here our model predictions become quite un- 
certain, as we lack sufficient information on the flux history of 
the years before the start of our campaign. There are only some 
sparse individual measurements available. 
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Fig. 14. Time evolution of the concentrations of several ions in components C and D. Further explanation see Fig. IT31 
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Table 6. Lower limits to the distance of the absorption compo- 
nents for different energy bands and instruments, at the 90% and 
99% confidence levels. 











90% lim. 


99% lim. 


Com- 


Instru- 


Range 


* 2 /d.o.f. 


logr 


logr 


ponent 


ment 






(m) 


(m) 


C 


pn 


0.5-1.5 keV 


508/279 


>18.26 


>17.99 


C 


pn 


0.62-0.78 keV 


109/59 


>17.84 


>17.62 


C 


RGS 


7-36 A 


7056/4972 


>18.29 


>18.11 


C 


RGS 


16.5-17.5 A 


217/169 


>18.32 


>18.11 


C 


RGS 


15.4-15.6 A 


45/19 


>17.58 


>17.37 


D 


RGS 


19.5-19.7 A 


33/19 


>17.16 


>17.02 


E 


pn 


0.5-1.5 keV 


508/279 


>17.15 


>17.01 


E 


RGS 


11.25-12.8 A 


342/259 


>16.50 


> 16.42 



5.4.2. Time evolution for the most prominent ions 

In Fig.[l4]we show the time evolution of the most important ions 
for our variability study (see Table |3j. Starting with oxygen, we 
see that the strong soft X-ray flux enhancement near the centre of 
the XMM-Newton observations has a strong ionising effect: both 
the O vii and O vm concentrations decrease for gas at distances < 
10 17 25 m; for longer distances, the history on longer time scales 
becomes more important. 

Comparing Fe x with O vn, it is seen that iron responds faster 
to continuum changes than oxygen: a similar delay is found for 
iron at three times greater distances (compare the curves for O vn 
and Fex at 10 1700 and 10 I7 5() m, respectively). For the iron M- 
shell ions Fe x - Fe xm the lower ionisations stages show a de- 
crease in column density while the higher ionisation stages show 
an increase near the centre of the XMM-Newton observations. 

Finally, the two panels for Fe xvn for components C and D 
show that the response not only strongly depends on the dis- 
tance to the central source, but also on the ionisation parameter: 
component D has a six times higher ionisation parameter than 
component C, and it responds in a less pronounced fashion. 

5.4.3. Changes in the transmission 

With the calculated time-dependent column densities, we next 
predict the changes in the transmission for each component as 
a function of distance. We compare the calculated transmission 
for individual observations to the calculated time-averaged trans- 
mission for all ten observations, similar to the spectral compar- 
isons shown in Fig. [6] This average is different for each adopted 
distance, cf. the different curves in Fig. [14] 

In Sect.|4]we found no evidence of an instantaneous response 
for any of the components where our sensitivity was sufficient 
to detect it. Here we follow a more general approach. For each 
ionisation component and for a grid of distances, we determine 
the formal agreement of the predicted ratio R pre d with the ob- 
served ,R b s by calculating ^ 2 -values in given energy or wave- 
length bands. We summarise our findings in Tabled In all cases, 
the best agreement was obtained for the assumption of no re- 
sponse (the low-density limit), so our data provide lower limits 
to the distances of the absorption components. The table also 
lists the ^ 2 -value and degrees of freedom for the low-density 
limit (no response). The offsets in ^ 2 -values from the degrees 
of freedom are caused by the slight mismatch of our continuum 
model for some of the observations (see Sect. 14.2b but these do 
not affect our distance estimates because they depend on the rel- 
ative increases in^ 2 and not on the absolute values. 



Table 7. Combined lower limits to the distance and density of 
the absorption components at the 90% confidence level. 



Component 


Distance 


Density 




(pc) 


(10" nT 3 ) 


C 


>71 


<0.28 


D 


>4.7 


<10.6 


E 


>4.6 


<1.7 



Table 8. Parameters for the model of a power law plus three 
Gaussians fitted to the pn spectra of Mrk 509 taken between 
2000-2006, divided by the 2009 spectrum as described in the 
text. 



yyyy-mm-dd f exp " Ar" Flux f N c d N D " 

2000- 10-25 20.4 0.296(8) 0.4244(16) -0.3(3) -1.3(4) -0.8(4) 

2001- 04-20 30.1 0.107(5) 0.6558(17) -1.8(4) -4.9(4) -2.5(4) 
2005-10-16 59.8 0.070(3) 0.7483(12) +0.6(3) +0.0(3) +0.0(3) 

2005- 10-18 32.4 0.006(5) 0.7885(17) -0.0(4) +0.2(4) -0.4(4) 

2006- 04-25 47.7 0.129(4) 0.7717(14) +0.1(3) -0.9(3) -0.6(3) 



Notes. Digits in parenthesis are the standard deviations in the last digits 
of the given value 

<n) Net exposure time (ks) of the pn observation; the total exposure time 
for the 2009 observations is 430.6 ks <6) Power-law index of the ratio; 
more positive numbers indicate harder spectra (c) Ratio of the normal- 
isations at 1 keV (</) Normalisations of the Gaussians for components 
C-E as labelled, multiplied by 100. 

For component C, the comparison of the pn results for the 
0.5-1.5 and 0.62-0.78 keV band shows that the highest sensi- 
tivity is reached by using the broader band. This is in line with 
Fig- El which shows that, apart from the main peak at 0.72 keV, 
there is also a weaker contribution from the broader 1 .0-1 .5 keV 
band. The RGS data in the 16.5-17.5 A band yield a distance 
limit that is almost equal to the limit obtained from the pn data. 
The narrower band around Fexvn (15.4-15.6 A) also yields 
lower limits, but these are less restrictive than the limits from 
the broader band. 

For component D, the pn data cannot distinguish the different 
cases well so we only get a lower limit from the O vm Lya line as 
measured by RGS. For component E the pn spectrum provides 
the best limit. 

Combining our results for component C from pn and RGS, 
and taking the most restrictive limits for components D (RGS) 
and E (pn), Table|7]summarises the 90% confidence lower limits 
to the distance and density of components C-E. 

6. Long-term variability 

In the previous section we have shown that there are no indi- 
cations of a change in the absorption components C-E during 
our intensive XMM-Newton monitoring period of 36 days. Here 
we investigate the longer-term variability. To this aim we have 
investigated the five archival XMM-Newton observations taken 
between 2000-2006. Similarly to the analysis of the data ob- 
tained during our campaign, we divided the fluxed pn spectra for 
all these observations by the fluxed pn spectrum of the combined 
2009 observations. As a first step, we simply fit the ratios of the 
fluxed spectra to a power law plus three Gaussians in the 0.5-1 .5 
band. For the Gaussians we keep the centroids and widths fixed 
to the values given in Table [2] representing the expected differ- 
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Fig. 15. Observed ratios of the fluxed pn spectra of Mrk 509 
taken between 2000 and 2006 to the stacked spectrum taken 
in 2009, in the 0.5-1.5 keV band. As described in the text, 
a local power-law approximation to the continuum ratios has 
been taken out. Dotted vertical lines indicate the centroids for 
three Gaussian components (Table|2| that represent the expected 
changes of the transmission associated to the ionising flux dif- 
ferences. The solid curves indicate the best-fit sum of three 
Gaussians to the residuals. For clarity of display, the curves for 
each spectrum have been shifted along the y-axis by multiples of 
0.1. 



ences due to changes in the ionising luminosity. We plot these 
ratios, with the local power-law model divided out, in Fig. [15] 
We list the best-fit amplitudes for the Gaussians in Table [8] 

It is clear from Fig. [15] and Table [8] that there are no sig- 
nificant differences in transmission between the 2005 and 2009 
spectra. There may be a hint for a small difference for compo- 
nent D in the 2006 spectrum, but it is not very strong. 

The strongest differences are visible in the 2000 and 2001 
spectra. For that reason, we have done a more sophisticated anal- 
ysis of the continuum spectrum for those two observations. We 
also extracted the XMM-Newton OM fluxes for those observa- 
tions and did spectral fits to both s pectra using exactly the same 
model as Mehdipo ur et al.l d201 ll) ; see also Sect. 15.21 for a de- 
scription of the spectral model they used. 

According to this continuum modelling, the ionising lumi- 
nosity during the 2000 and 2001 observations is lower by -0.45 
and -0.23 dex. We calculated the corresponding difference in 
transmission using the method described in Sect.|4]for the 2009 
spectra. We show the results as a comparison of the observed 
ratio /? b s to the predicted ratio R pie d (cf. Eq.[2]-[3]l in Fig. [TBI 

The broad dips that are visible in the 2000 and 2001 spec- 
tral ratios shown in Fig. [15] between 0.7-1.2 keV have almost 




1 1.2 
Energy (keV) 



2001-04-20 




Energy (keV) 

Fig. 16. Data points with error bars: observed ratios of the spec- 
tra taken in 2000 (upper panel) and 2001 (lower panel) to the 
average spectrum taken in 2009 in the representation of Eq. (0. 
We use here realistic continuum models (see text) rather than 
simplified local power laws as used in Fig. [15] to take out the 
changes of the continuum. Solid green lines: predicted ratio in 
the representation of Eq. (O for ionisation component C; dashed 
blue lines: same for component D; dotted red lines: same for 
component E; thick magenta line: components A-E combined. 
The predictions for components A and B are close to zero and 
not shown here separately. 



disappeared. This is due to the more realistic broad-band optical 
to X-ray modelling of the continua in the present case, com- 
pared with the local 0.5-1 .5 keV power-law approximation used 
in Fig. [15] Furthermore, the spectra are not sensitive enough to 
detect changes for component C. For component D, both spectra 
show a narrow dip in the data around 0.8 keV. When fitted with a 
Gaussian with centroid and width frozen to the values for com- 
ponent D taken from Table [2] we find best-fit amplitudes TV of 
-0.014 + 0.006 and -0.020 + 0.004 for the 2000 and 2001 spec- 
tra, respectively. These values correspond to a nominal decrease 
in <f of -0.23 + 0.10 and -0.33 + 0.07 dex. Finally, it is evident 
that the large changes that are predicted to occur for component 
E are not present. Thus, component E appears to be stable in ion- 
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isation parameter over ten years. This can only be explained by 
a low density, such that either there are no changes at all over ten 
years, or that the average ionising fluxes over the few months be- 
fore each of the 2000, 2001 and 2009 observations were similar 
by coincidence. 

Comparison of the 2000 and 2001 spectra shows that the 
quality of the data does not allow us to establish whether there 
are differences in the transmission between these spectra. 

7. Discussion and conclusions 

In this paper we have constrained the distance of the dominant 
X-ray outflow components in Mrk 509. This has been a chal- 
lenging task, given the complexity of the problem and the refined 
analysis that is necessary. 

The key to achieving these results was the deep 600 ks time- 
averaged RGS spectrum that allowed us to derive the ionisa- 
tion s tructure of the outflo w as seen through X-ray absorption 
lines (iDetmers et al.ll20TTl) . The discrete components found in 
that work allowed us to make predictions for the time-dependent 
transmission of the outflow for different distances to the cen- 
tral source, and the distance limits were obtained by comparing 
those predictions with the observed spectra. 

Another necessary ingredient in our analysis was the accu- 
rate broad-b and spectral modelling from the optical to the hard 
X-ray band (Mehdinour et al.11201 ll) . This was not only impor- 
tant for deriving the shape of the ionising continuum, but also 
for obtaining the underlying continuum spectrum for each ob- 
servation in the energy band around 1 keV where most of the 
changes in transmission from the outflow are expected to occur. 

We have used here the continuum models as derived by 
iMehdipour et al.l (1201 ll) . These were based on the measured OM 
and pn spectra, but also on the scaled, non-simultaneous FUSE 
spectra described in that paper. This gives an excellent empirical 
description of the spectra (see also below). Relieving the con- 
straints from FUSE yields somewhat different parameters for the 
continuum components from the ones used here (see Petrucci et 
al. 2012, in preparation), but their effects on the present results 
are minor. 

Our analysis has shown that in the pn-band the continuum 
model spectra are accurate down to the 1% level (Sect. 14. 2t 
and that narrow-band time variations can be detected with the 
pn down to the 0.5% level (Sect. 12.1.11) . For our Mrk 509 data, 
an instantaneous response to the ionising continuum would pre- 
dict changes up to 2% at certain energies for some of the ionisa- 
tion components (for example between spectra 1 and 4, Fig. [6]). 
Given the detection limits that we mentioned, this is certainly a 
detectable signal. For the somewhat lower ionised components 
with the dominant absorption troughs below ~ 1 keV, RGS gives 
useful constraints through the same blends as visible with pn, 
but also through some strong individual absorption lines (Fe xvn, 
O vm). 

However, our data show no evidence of an instantaneous re- 
sponse to the continuum variations, as outlined in Sects. [4721 and 
14.31 Therefore we exclude the high-density limit for the gas, im- 
plying that the absorber must be located farther away from the 
central black hole. Our monitoring with S wift before the onset 
of the actual XMM-Newton monitoring dKaastraetalJl2011bt 
Mehdipou r et al.ll201 lb has been very important. It allowed us 
to develop and use a realistic time-dependent photoionisation 
model that predicts the change in transmission for a broad range 
of recombination time scales, hence distances to the ionising 
source. 



In our time-dependent calculations, we assumed equilibrium 
for the first Swift observation as an initial condition. Because 
the first XMM-Newton observation starts six weeks later, most 
uncertainties due to this initial condition will be washed out for 
time scales shorter than this. On longer time scales this is not the 
case. We note, however, that the flux during our XMM-Newton 
observations was slighty higher than the long-term average flux, 
while the first Swift observations was perhaps closer to this long- 
term average. 

Again comparing our predictions with the data, we have put 
tight lower limits to the distance and density of components C- 
E. Comparing our spectra taken in 2009 with archival data, in 
particular those taken in a relatively low flux state in 2000 and 
2001, showed changes for component D on time scales of 8.5 
years. This is 50 times the duration of our monitoring campaign, 
so it yields a 50 times lower limit to the density compared to the 
number given in Table|7] of 2.1 x 10 8 m -3 , with a corresponding 
upper limit to the distance of 33 pc. 

In the long term, we see no detectable changes for com- 
ponents C and E. This does not imply automatically that no 
changes have occurred for those components, but we lack ade- 
quate monitoring before the 2000 and 2001 observations. Given 
the lack of short-term variability for components C-E, the ion- 
isation equilibrium depends on the time-averaged spectrum for 
periods of months to years. 

Because Mrk 509 can show ionising flux variations up to 
about 40%, this means in practice that for any spectrum taken, 
the ionising flux used to calculate the ionisation balance for the 
analysis of the absorber is also uncertain by that order of mag- 
nitude. In our case, the average flux during our XMM-Newton 
campaign may have been somewhat enhanced, perhaps by 15- 
20%, compared to the long-term average, due to the outburst 
that peaked around the epoch of spectra 4 and 5. This would, 
in retrospect, justify some minor corrections to the precise val- 
ues of $ that we der ived from the time-average RGS spectrum 
(IDetmers et al.ll201lT) . but they do not alter the general picture. 

It is possible to derive upper limits to the distance using a 
rather simple argument. Our lower limits to the distance r im- 
ply an upper limit to the density n through n = L/^r 2 . Because 
the typical thickness of the absorbing layer a is related to the 
measured column density A^h = na, we derive a lower limit to 
the relative thickness a/r. Reversing the argument, because a/r 
must be smaller than unity, we get an upper limit to the distance. 
For our components C-E these are 9 kpc, 1.3 kpc, and 21 pc, 
respectively. The limit for component C is not very useful, and 
for component D we have already found an upper limit of 33 pc 
from the long-term variability, but the limit for component E is 
interesting. 

There is a caveat here, however. In this paper, we have used 
the five ionisation components A-E der ived from a d i rect fit 
to the total measured column densities. Detmers et al] (1201 ll) 
also fitted the RGS spectrum directly using the xabs model of 
SPEX. The components found that way are, within their sta- 
tistical uncertainties, compatible with what we use here, but 
for component E there are larger uncertainties. We have used 
log£ = 3.60 + 0.27, while the direct fit with t he xabs model 
yielde d log£ = 3.26 + 0.06. As discussed by Detmer set al.l 
(2011) this is the same component; the differences arise mainly 
because component E is dominated in the RGS band by Fe xxi 
and Fe xxn. Making the ionisation parameter larger but also in- 
creasing the total hydrogen column yields almost the same ionic 
column densities for these ions. 

We re-ran all our procedures but with the ionisation param- 
eters and column densities replaced by those based on a direct 
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xabs fit (Table 5 in Detmers e t al.l201 H) . In particular for compo- 
nent E the upper limit to the distance from the constraint ajr < 1 
increases significantly, from 21 pc to 400 pc. For the other com- 
ponents the differences are not very important. 

Using the Ch andra LETGS data of our campaign, 
lEbrero et alj d201 lb also derived upper limits of 300 and 50 pc 
for the components C and the blend of D and E, respec- 
tively. These limits are consistent with our present findings, but 
Ebre ro et alj d201 lb had to make the additional assumptions that 
the solid angle of the outflow is ~ ji sr, and that the momentum 
of the outflow is clo se to the momentum of the scattered and ab- 
sorbed radiation, cf. lBlustin et al.1 d2005l) . Our present estimates 
are more direct. 

iDetmers et alj d2010l) analysed the XMM-Newton spectra of 
Mrk 509 taken in 2005 and 2006. A simultaneous fit to the pn 
and RGS data of these spectra showed three ionisation compo- 
nents that correspond to our components B, C, and a blend of 
D and E. Components B and C did not show any significant 
variab ility between the observations. The fits of IDetmers et alJ 
(2010) indicate that the ionisation parameter for component D/E 
had decreased from 3.29+0.04 in the first observation of 2005 to 
3.01+0.06 in the 2006 observation. From this change within 0.5 
year an upper limit to the distance of 0.5 pc was derived. Such a 
value clearly contradicts our lower limit of about 5 pc for compo- 
nents D and E. With our present model, however, a change in 0.5 
year would correspond to a distance of about 15 pc, much larger 
than 0.5 pc. The discrepancy can be traced dow n to the use of 
Fe xn t o estimate the recombination time scale by Detmer s~et al.l 
d2010l) . This ion has its peak concentration near the ionisation 
parameter for component D/E. It is therefore relatively insen- 
sitive to changes in the ionising flux. As a consequence, the 
formal recombination time scale (Eq. 4 in Det mers et alj feOlO) 
becomes very large, leading to high density and short distance. 
Furthermore, the apparent change in transmission between 2005 
and 2006 is possibly affected by some systematic effects and 
less significant than the formal 3.9<x significance quoted above. 
Factors contributing to a larger uncertainty are the lack of multi- 
pointing mode data for the RGS spectra in 2005 and 2006, the 
shorter exposure times leading to the detection of only three out 
of the five ionisation components, blending of ionisation com- 
ponents D and E, and global fitting with imperfect effective area 
instead of using ratios like in the present paper. The combined 
difference signal from components D and E between the first 
2005 and the 2006 spectrum taken from Table [8] are, however, 
still consistent with a change at the 2.5cr significance level. We 
keep it open here whether this change is real or not, but if it is 
real, the upper limit of ~15 pc, according to our improved cal- 
culations, is close to the limit of 21-400 pc obtained from the 
requirement that a/r < 1 for component E as derived earlier in 
this section. 

The distance limits derived in this paper (Table |7]i imply an 
origin in the NLR or torus region of Mrk 509, rather than in the 
accretion disc. The limits are fully compatible with the upper 
limit of 250 pc deri ved from UV va riability of the fastest UV 
outflow component dKriss et al.ll201 ll) . This does not exclude an 
accreti on-disc wind com ponent. Perhaps the highly ionised, fast 
(0.05c iPonti et al.ll2009l). o r ultra-fast (0.14-0.20c. lCappi et al l 
120091: iTombesi et al.1 1201 Oh outflows found in archival high- 
energy spectra of Mrk 509, but not found in our present data, 
could be examples of such a disc outflow. We defer a comprehen- 
sive overview of the various absorption components, combining 
both the ultraviolet and X-ray data, to another paper (Ebrero et 
al. 201 1, in preparation). 



In summary, we derived distance limits for the main X-ray 
absorption components in Mrk 509. Component C has a distance 
of > 70 pc, component D is between 5-33 pc, and component E 
has a distance > 5 pc but less than 21^00 pc, depending upon 
modelling details. 
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